Instrinsic oscillations of treadmilling microtubules in a motor bath 
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We analyse the dynamics of overlapping antiparallel treadmilling microtubules in the presence 
of crosslinking processive motor proteins that counterbalance an external force. We show that 
coupling the force-dependent velocity of motors and the kinetics of motor exchange with a bath 
in the presence of treadmilling leads generically to oscillatory behavior. In addition we show that 
coupling the polymerization kinetics to the external force through the kinetics of the crosslinking 
motors can stabilize the oscillatory instability into finite-amplitude nonlinear oscillations and may 
lead to other scenarios, including bistability. 
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The collective action of motor proteins plays a funda- 
mental role in biological self-organization [l|, . Sponta- 
neous mechanical oscillations in biological systems where 
groups of motors are elastically coupled to their environ- 
ment constitute a generic phenomenon Q. Examples of 
such oscillations can be found, e.g. in the coupling of 
the mitotic spindle to the cell environment [1, or the 
mechanical oscillations in muscle fibers [gl; the oscilla- 
tory behavior of motor-filament complexes has also been 
invoked to play a crucial role in mitosis Two main 
simplifications are usually considered in most studies of 
collective action of motors coupled to filaments, such as 
in gliding assays or other in vitro experiments [3[ : (i) the 
treadmilling dynamics of the filaments is eliminated; (ii) 
motors are fixed onto a substrate. In biological situa- 
tions, however, it is common that filaments are contin- 
uously treadmilling and motors can freely be exchanged 
with a motor bath in solution A paradigmatic 

example are the interpolar microtubules of the mitotic 
spindle, which undergo the so-called poleward flux [gj. 

In this letter we show that force-dependent velocity 
of processive motors crosslinking treadmilling filaments 
against an external constant force, leads generically to 
intrinsic oscillations, which may be either damped or am- 
plified. We also show that, in the latter case, the oscil- 
latory instability can be stabilized into nonlinear oscilla- 
tions by the coupling of kinetics of motor-exchange with 
microtubule polymerization. 

To elucidate the generic character of the phenomenon, 
we focus on a minimal arrangement, inspired by the inter- 
polar microtubules in the mitotic spindle, consisting of a 
symmetric pair of polymerizing antiparallel microtubules 
(MTs). These are crosslinked by (processive) tetrameric 
kinesins (Eg5) in a central region of overlap of size i, 
typically much smaller than the MT length The 
MTs overlap in their plus ends, the extremum where they 
polymerize as they are being pulled by the motors. By 
symmetry, the crosslinking motors remain in the overlap 



region while they slide the two MTs in opposite direc- 
tions towards their respective minus ends (the so-called 
poleward flux of spindle MTs). At the same time, the 
crosslinking motors balance an external inward force F 
that is assumed constant. 

In Ref. [lo| it was shown that under some plausible 
assumptions, the dynamics of the overlap region can be 
decoupled from that of the MTs outside, which enters 
through boundary conditions. Then, the motor kinet- 
ics in the overlap region may be described by simple 
rate equations for two populations of motors in the over- 
lap region, ric designating the average number of motors 
crosslinked to both MTs, and ni, those bound to either 
one of them. The kinetic equations take the simple form 

= -2ku{nc)nc + hrit (1) 

at 

^ = 2kMn,-{ki, + Oni, + kl''p3Di, (2) 

where kf, and k^ are respectively the binding and un- 
binding rates of the motor domains for motors that are 
not under load. Only crosslinked motors exert forces 
and we assume a force-dependent kinetics of the stan- 
dard form fc„ = k^exp{fmb/kBT), where 6 is a length 
in the nanometer scale characterizing the activation pro- 
cess, ksT is the thermal energy, and the force per mo- 
tor is taken as fm = F/ric. This results in a strongly 
nonlinear dependence of k^ on ric. We keep this force- 
dependent kinetics to remain quantitatively as realistic 
as possible although this dependence is not necessary to 
explain the qualitative picture of emergence of oscillatory 
behavior. The exchange with the motor bath is described 
by the balance between the last two terms in Eq. ([2]) JJl . 
The motor intake from the bath is controlled by the rate 
k^^ and is proportional to the motor density in solution 
PsD, which together with F deflne the two experimen- 
tal control parameters of the problem. In general, the 
overlap length, will evolve dynamically as a result of 



2 



the balance between the overah polymerization, V^, and 
sliding, Vs, velocities 



d£ f F 

— = 2{Vp - Vs) = 2FpK) - 2Vb 1 y 
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(3) 



thus relaxing a strong assumption of Ref . [l^l ■ For sim- 
plicity and consistently with experimental evidence, we 
have assumed a linear decrease of motor velocity relative 
to the MT with the force per motor /„ — F/ric, and 
where is the single-motor stall force. The dependence 
oiVp{nc) expresses the coupling between motor dynamics 
and polymerization. Since the emergence of the oscilla- 
tory behaviour itself does not depend on this coupling, 
for the first part of the analysis we will take Vp = const. 
Later on we will discuss how this additional coupling may 
actually stabilize nonlinear oscillations. 

Eqs. ([J)-© provides the essential description of the 
motor-microtubule complex. They can be expressed in 
dimensionless form in terms of the characteristic detach- 
ment force, kBT/b, the processivity length, Vb/fc°, and 
the detachment rate, fc^. Accordingly, we will use the 
dimensionless variables F — Fb/KsT, f = fsb/KsT, 
he = njF, hb = Ub/F, £ = ifk°/Vo and r ^ tk^. We 

define A = ^^^ttt^tt^j a ratio of the strength of motor 
influx and outflux in the overlap region, as a dimension- 
less parameter that expresses the balance between the 
two external control parameters that couple the system 
with the environment. In turn, 7 = kb/k^ measures the 
asymmetry in motor attachment/detachment at vanish- 
ing load. The relevant parameters become then A, 7 and 
g = /(I — Ip/Vo). Eqs. (HI)-© have a stationary solution 
of the form 
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J, - — exp[g], = —r exp[g] - 2g, (4) 
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which holds for g > {Vp < Vq). Physically meaningful 
configurations are further restricted to the regime A < 
^ to ensure ^ > 0. 

The linear stability of the fixed point is determined by 
the eigenvalue equation 

+ (1 + 7 + G)A2 + {2g^A - G)X + 2g^A = 0. (5) 

where we introduce G = (g — 1) exp{g) for simplicity. 
The conditions for the existence of three real roots, and 
the corresponding regions of stability can be found ex- 
plicitly, and will not be discussed here. We focus on 
the more interesting case where two complex conjugate 
modes appear, 



Ai = -1 -7 + G 
A2,3 = Ao ± i9, 



(6) 
(7) 



where Aq is real and is known explicitly in terms of 
the model parameters, Ai is always negative, and 6 = 



y^2g^A — G. The stability region corresponds to Aq < 0, 
while the stability boundary Ao = is defined by 
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2g'^A-G' 



(8) 



where stability also requires 2g'^A — G > 0, imposing a 
lower bound on A. This constraint, together with ^ > 
does not allow for stable, physical steady state configura- 
tions when g > 2. Since 7 and A are positive parameters, 
the steady antiparallel arrays are intrinsically stable for 
g < 1, regardless of A. Only in the regime 1 < g < 2 the 
motor/MT complex stability is controlled by A. Fig. [1] 
displays the stability curve, as a dashed line, for g = 7/4 
where stable antiparallel structures exist in a bound re- 
gion of parameters, Am < A < Am- The instability 
associated to the change of sign of Aq corresponds to a 
Hopf bifurcation, giving rise to damped or growing os- 
cillations at a finite frequency respectively at the stable 
and unstable sides of the stability boundary. The char- 
acteristic dimensionless frequency 9 = Im{\2) is given 
in full units by ujb = k^\j2g'^A — G and depends both 
on kinetic parameters and the environment through A; 
the oscillatory behavior persists outside the neighbour- 
hood of the instability, with a weak dependence of the 
frequency on nonlinearities. 

To gain further analytical insight, it is useful to focus 
on the region near the instability threshold, where | Ai | 3> 
|Ao|. We can then write Aq — G — 7 and Ai ~ — 1-l-Ao and 
exploit the separation of time scales between the (slow) 
oscillatory modes and the (fast) relaxation mode. This 
is a rigorous two-dimensional reduction of the problem 
(to the so-called center manifold) which can be extended 
to the nonlinear level. The adiabatic elimination of the 
fast mode 5(t>i ^ Gdhc + j6hb — 9^6£/2g^, which decays 
as S(f)i ^ exp(— Aii), implies that it will be slaved to 
the oscillatory modes. It is therefore legitimate to take 
6(j)i ~ which introduces a constraint on the variations 
of the three variables. Hence, at linear level there is an 
approximately conserved quantity of the form 



1 9"^ 
j{Gnc + inb) ^ 



(9) 



This implies that the overlap length, reacts to keep 
constant an effective motor density that weights the 
two types of motors differently; only when A ^ 1 this 
weighted density reduces to the total motor density. Con- 
sequently, as long as |Ai| » |Ao| ~ |G — 7I, the dynamics 
of Sric at linear level is that of a harmonic oscillator 



6fic — \5hc + 9 Sh 



0. 



(10) 



where the oscillations of frequency 9 are exponentially 
damped or amplified depending on the sign of Ao. The 
other two variables have the same behavior with differ- 
ent amplitudes and phases which follow from the linear 

relation d£ = —2g^Shc. 
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FIG. 1: Stability diagram of MT array: The stability line 
for a* = (dashed curve) is derived from Eq.([8}, where no 
limit cycle behaviour is observed. For a* — 0, the MT array 
is unstable below Am = 'oe'^'^ jT^ and above Aa/ = 6^/2"^ 
(as MT arrays are unphysical with / < 0). The area en- 
closed within the dashed curve is the region of stable arrays 
for a* = 0. Stability region is enhanced for inhibitory motors 
with a* = 0.01. The smooth curve (red online) depicts the 
stability line for this case. Between the smooth and dotted 
curves (maroon online), the arrays are nonlinearly stabilized 
and limit cycles are observed. Here g = \- 



It is important to remark that the overlap length I 
must be coupled to the motor dynamics in two ways in 
order to produce oscillations: (i) through the motor ex- 
change in the bath with a motor intake proportional to 
^; and (ii) through the decreasing velocity-force curve of 
the motors, implying that ^ —Sric- The existence of 
oscillatory behavior is independent of the detailed depen- 
dence of motor kinetics and motor velocity on the applied 
force (as long as it is monotonously decreasing). The 
kinetic exchange between and Uc is nevertheless es- 
sential to control the existence of a nontrivial stationary 
point (with finite values of all variables) and its stability, 
in particular allowing for unstable growth of the oscilla- 
tions. A nonlinear analysis is necessary in that case, and 
model details may come into play. 

For constant Vp, the amplitude of the oscillations is not 
saturated by nonlinearities leading eventually to array 
disassembly. The coupling to the external environment 
or additional dissipation provides a means to stabilize 
the filament-motor assemblies. Remarkably, simple as- 
sumptions on the coupling between motor kinetics and 
MT polymerization provide possible mechanisms of non- 
linear stabilization of the oscillations in some regimes, in- 
troducing new dynamical scenarios. As an illustration we 
may consider Vp{nc) = [Cp{nc) — Cd{nc)] e where Cp and 
Cd are the polymerization and depolymerization rates 
at the plus end, with e the length increase per added 
monomer. If the coupling is weak, one can linearize 



the dependence of motor concentration in the polymer- 
ization rates, rendering the detailed functional form of 
the coupling irrelevant. Accordingly, one can generically 
write Cp((i) ~ C'p(dp)(l ~ Q!j5(dp)^c/^)- In this weak cou- 
pling limit, it is enough to consider a single parameter, 
ao, the coupling strength, which characterizes the rel- 
ative magnitude of the polymerization and depolymer- 
ization rates motors promote depolymerization for 
ao < and inhibit it for ao > 0, and their effects can 
be quantified in terms of the dimensionless parameters 
a = aoFf/lp, Ti = eCp/ik°lp) and T2 = eCdpl[k%), 
giving now .g = /[I — (Fi — r2)]. Eq. ([3]), takes now the 
dimensionless form 
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where a* — fT2a. The steady solution then reads, 
g ~ l/ric A Uc Uc 



(11) 



(12) 



along with fih = ^exp(J-). For a* ^ ^^.9^, the steady 

antiparallel array changes gradually with fi^ = l/.g + 
a*/€/g3_ 

The nonlinear dependence between Uc and £ allows for 
physical solutions if nc,m > I/5, £m < ^01* / g'^ and 1/A > 
a/g, which imply g > 2. As a result, qualitatively new 
steady arrangements appear, unaccessible for a = . For 
g > 2 two new steady configurations are allowed, and the 
overlap region has always a positive length regardless of 
A. The numerical analysis shows that polymerization 
enhancement, a* > 0, favors the array stability, as shown 
in Fig. [T] for a* — 1/100, while the opposite holds if 
motors enhance MT depolymerization jl5| . 

For inhibitory couplings, a > 0, linearly unstable ar- 
rays can be stabilized nonlinearly close to the stabil- 
ity curve leading to limit cycle oscillations. The ampli- 
tude of these nonlinear oscillations increases when mov- 
ing into the unstable region, leading eventually to ar- 
ray disassembly. For a — 1/100, Fig. [T] shows the re- 
gion where nonlinear oscillations are sustained. Fig. ^a.) 
identifies a regime, A > 0.8 where the amplitude van- 
ishes as the square root of the distance to the linear 
stability threshold, analogous to a second order transi- 
tion while for A < 0.8, the nonlinear stabilization has 
a finite amplitude from the outset, reminiscent of a first 
order transition, as shown in Fig.[2jb).The oscillation fre- 
quency always varies as we move away from the instabil- 
ity threshold, as depicted in Fig. [2j As one approaches 
the nonlinear stability threshold, decreasing 7, the fre- 
quencies decrease. We have never observed nonlinearly 
stabilized arrays when motors promote MT depolymer- 
ization (a < 0). 

The existence of three steady configurations for g > 2 
when a > allows for new dynamic scenarios. As shown 
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FIG. 2: (a) Variation of Ai, the amplitude of tiie non- linear 
oscillations of the overlap region, I, with 7 for a* = 0.01, g = 
1.75. Circles (red curve online) corresponds to A = 0.9, when 
the instability is supercritical and Squares (maroon curve on- 
line) for A — 0.74, when it is subcritical. b) Variation of 
the frequency of the limit cycles as a function of 7 for the 
corresponding points. 
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FIG. 3: Stability diagram of an antiparallel MT array for g — 
2.7 and a* = 0.1. At large 7 there is a region of intermediate 
A values where both a short (1) and a long (2) MT arrays are 
stable, leading to bistability. 



in Fig. [21 by varying A it is possible to find bistabil- 
ity. The coupling between the array nonlinear oscilla- 
tions with different external frequencies in this bistable 
regime provides new scenarios to control the motor-MT 
complexes stability and its sensitivity to changes in the 
array environment. 

We have shown that the self-organized, coupled dy- 
namics of the overlapping region of antiparallel MT ar- 
rays gives rise to intrinsic oscillations due to the local- 
ized motor kinetics in the overlap region and its coupling 
with MT polymerization. These oscillations may be rel- 
evant in spindle dynamics, where oscillations reported to 



date emerge from the interaction of the antiparallel array 
with the surrounding media through centrosomal micro- 
tubules or the interaction of molecular motors with 
the chromosomes MT attach to 0. The reported oscil- 
lations are generated at relatively small scales, of the or- 
der of the few fractions of a micrometer which character- 
ize the overlap between biofilaments, but the frequency 
range, controlled by the motor unbinding rates, is typ- 
ically of the order of s~^. The coexistence of different 
stable arrays also allows for bistability between antipar- 
allel arrays, providing enhanced mechanical versatility of 
these structures to environmental changes. 
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In Eq.((TJ a term of incoming bound motors from the 
non-overlaping region such as that considered in Ref. 
could be introduced, given the processivity of motors. For 
simplicity, and since such effect is relatively small we omit 
this term here. 

In-vitro experimental evidence shows that processive 
plus-end motors enhance MT depolymerization in bud- 
ding yeasts pTi . [l^ . Nonetheless, the interactions of 
molecular motors with MT in other situations, and the 
possibility that other molecules may modify the effective 
motor/MT interactions leaves this issue generically open. 
In our approach all these details are effectively included 
in the phenomenological parameter a. 
When motors promote MT polymerization the restriction 
g < 2 holds. In this case A^ and Am approach as a de- 
creases, reducing the parameter range where antiparallel 
arrays are stable. 



